DO-TH-02/17 

October 2002 
revised January 2003 



The 0(N) linear sigma model at finite temperature 
beyond the Hartree approximation 

Jiirgen Baacke 1 and Stefan Michalski 2 

Institut fur Physik, Universitdt Dortmund 
D-^221 Dortmund, Germany 



Abstract 

We study the 0(N) linear sigma model with spontaneous symmetry breaking, using 
a Hartree-like ansatz with a classical field and variational masses. We go beyond the 
Hartree approximation by including the twodoop contribution, the sunset diagram, 
using the 2PPI expansion. We have computed numerically the effective potential at 
finite temperature. We find a phase transition of second order, while it is first order 
in the onedoop Hartree approximation. We also discuss some implications of the fact 
that in this order, the decay of the sigma into two pions affects the thermal diagrams. 
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1 Introduction 



The O(N) linear sigma model has a long-standing history, in particular as a basic model 
for a quantum field theory with spontaneous symmetry breaking PJ El El Ej Early inves- 
tigations beyond the classical level have been based on including one-loop quantum and 
thermal corrections. These studies have been centered around the discussion of the one- 
loop effective potential Kfr(0) where <fi is the mean value of the quantum field in a sense 
being defined more precisely by the effective action formalism, summing up one-particle 
irreducible (1PI) graphs. A next class of approximations include bubble resummations, as 
motivated by the large- JV limit. In the model with spontaneous symmetry breaking one 
finds a second order phase transition such that the symmetry is restored at high tempera- 
ture. 

Another approximation, going somewhat beyond the leading order of the large- JV ex- 
pansion, is the Hartree approximation; it includes only local one-loop corrections to the 
effective mass and thereby takes into account some, but not all, next-to-leading order cor- 
rections in 1/JV. The Hartree approximation of the O(N) linear sigma model has been 
studied at finite temperature by various authors El 13 El El HH ITU] : the model with 
spontaneous symmetry breaking is found to have a phase transition of first order towards 
the symmetric phase at high temperature. In contrast to the large- N case the mass of the 
pion quantum fluctuations does not vanish in the broken phase. This has been discussed 
as a "violation of the Nambu-Goldstone theorem" (see e.g. Ref. jH] and references therein); 
the presently accepted point of view [HI El E] is that the "sigma and pion masses" in the 
Hartree scheme are just variational parameters, and not the real pion and sigma masses, 
which are to be computed from the effective potential at its minimum. 

If one wants to go beyond the large- N and Hartree approximations there is a variety 
of choices. The systematic expansions are based on the resummation scheme by Cornwall, 
Jackiw and Tomboulis (CJT) the 2PI scheme. Within this scheme one may select 
certain groups of graphs in order to obtain systematic expansions in 1/JV or in the number 
of loops (order in Ti). Beyond the leading order, these extensions require technically quite 
involved analytical and numerical calculations [12 j . In general one has to solve Schwinger- 
Dyson equations for the Green functions which in the present case would even form a 
coupled system. Little is known about the merits of the next-to-leading order extensions 
as such calculations at finite temperature in 3 + 1 dimensions are not yet available. 

A technically less demanding approach is the 2PPI resummation introduced by Ver- 
schelde |14[ll5j. Here, instead of treating the Green functions as variational parameters one 
just introduces variational masses, like in the Hartree approximation. This implies that 
the resummation is only over local insertions, the 2-particle "point reducible" graphs, i.e., 
graphs that fall apart if one cuts two lines meeting at the same point (the 2PPR point). 
This approach is based on a variational principle for expectation values of local composite 
operators, i.e., all of the system's equations of motion can be derived from a single func- 
tional. The 2PPI effective action is identical to that in the Hartree approximation if only 
one-loop 2PPI graphs are included; this has been studied in Ref. [TH ITU] . For the complete 
two-loop approximation one must also include the sunset diagram. For the case N — 1 
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Smet et al. ^H] have evaluated the effective potential; they found that instead of a first 
order phase transition one obtains a second order one. Here we extend this investigation 
to the case of general N. 

The O(N) linear sigma model with spontaneous symmetry breaking has been studied in 
nonequilibrium quantum field theory as well, mostly in the large- N limit and with different 
initial conditions for the mean field <fi = ($) and for the density matrix of the fluctuations 

[mim[3l23inil22l23l211l2S]- There the mean field <f> = ( $ ) becomes time dependent. 
As far as symmetry restoration is concerned striking similarities with finite temperature 
quantum field theory are observed [El [TBI ES] : If the system is supplied with a high 
initial energy density it displays symmetry restoration at late times in the sense that the 
mean field settles at <fi = or oscillates around this value, while at lower energy densities 
the system ends up in a broken symmetry phase where the time average of (p(t) remains 
different from zero, and where the pion mass, the time-dependent mass of the quantum 
fluctuations, goes to zero. This phase structure persists if one uses the Hartree instead of 
the large- N approximation [26J; in this case, as in thermal equilibrium, the effective mass 
of the pion fluctuations remains finite even at low energy densities. 

However, in the large-iV or Hartree approximations the system does not approach 
thermal equilibrium. This problem has been addressed in a general way in Ref. [2*7[ 
EEj. Numerous authors [23 [SOI EH IH21 EHl EH EHl EE] have tried recently to find useful 
approximations beyond the leading orders. Up to now numerical simulations are mostly 
limited to 1 + 1 dimensional models. Most of the new approximations show large deviations 
from the large- N approximation, and they indicate thermalization. The proper case of an 
O(N) model with spontaneous symmetry breaking has not been investigated up to now. 
Indeed higher corrections have not even been included in equilibrium calculations for such 
models. If one tries to appreciate the quality of various approximations such equilibrium 
computations should be able to yield useful additional insights. It is one of the purposes 
of this work to initiate such investigations. 

The plan of the paper is as follows: in section [21 we present the general formulation 
of the model and of the 2PPI formalism. In section El we explicitly formulate a potential 
U(m^,m^, 0) that by variation of and leads to the gap equations. The technical 
details of the relevant Feynman graphs and a comparison of the 2PPI expansion to CJT's 
2PI approach are presented in the appendices. In sectionHJwe discuss our numerical results, 
we end with a summary and an outlook in section El 

2 Basic equations 

The Lagrange density of the O(N) linear sigma model is given by 



where <3>j is a vector with N components. We intend to compute the effective potential of 
this model at finite temperature. This model has been studied at large N and in the Hartree 
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approximation, which both represent bubble resummations. One of the possibilities to go 
beyond these approximations, and in particular to include higher loop corrections is the 
use of the 2PI or CJT formalism; this is technically involved, even in equilibrium, as one 
has to solve Schwinger-Dyson equations for the Green functions, in the present case indeed 
a coupled system of integral equations. 

Another possibility of going beyond the leading order approximations has been proposed 
by Verschelde |14[ I15j. the so-called 2PPI formalism. This is a variant of the 2PI (CJT) 
formalism by Cornwall, Jackiw and Tomboulis [T^]. In the 2PPI approach the composite 
operator is local while in 2PI it is bilocal. Here the resummation encompasses all 
2-particle point reducible graphs, graphs that fall apart if two lines meeting at one point 
(vertex), the 2PPR point, are cut. These graphs are deleted in the 1PI effective action, 
which thereby is replaced by the 2PPI effective action. They are taken into account by a 
mass insertion like in the Hartree approximation — to which the 2PPI expansion reduces 
in the one-loop approximation. We compare this approach to the well-known 2PI CJT 
formalism in Appendix [X] 

The problem occurring in the Hartree approximation, namely the lack of a consistent 
renormalization, has been solved in a systematic way. The inconsistencies are avoided by 
recognizing that in the resummation the counterterms have to be divided into 2PPI and 
2PPR parts. The 2 particle point reducible parts renormalize the gap equation, the 2PPI 
parts renormalize the 2PPI effective action. This procedure has been discussed in technical 
detail in [THl HHj and been applied to a first two-loop calculation for the N = 1 model (THj , 
including the sunset diagram as the only 2-loop 2PPI term, the only other new terms being 
one-loop graphs computed with the one-loop counterterm Lagrangian. 

We will not go into details here. For the O(N) case we use the explicit formulae of 
[T5] . The classical field is denoted by (pi = the bubble resummation is defined by 

introducing local insertions Ay = — 4>i4>j = {&i&j)conn. which collect all 2PPR 

graphs. The resummation is defined by including these insertions as well as the seagull 
insertions, it is obtained by introducing into the propagators the effective mass 3 

m% = -Xv^ij + 2A [(Mj + A -) + A (0 fc fc + A* fc )] . (2.2) 

The motivation and formal derivation of the 2PPI effective action for the case N = 1 is 
presented in Appendix El The generalization to arbitrary iV is straightforward ^U]. The 
2PPI effective action can be written as 

T = S dass + Tf p %, mj] - - (AjjAjj + 2AyAii) . (2.3) 

It includes all 2PPI graphs as defined above, with the mass terms replaced by the variational 
masses m^, and it is computed using the 2PPI parts of the counterterms. The last term 
is introduced in order to avoid double counting. The local self-energies A^- can be shown 
to be related to the "quantum" part of the 2PPI action via 

1 dT 2PPI (m 2 -) 

\^ = q jy (2.4) 

2 dmjj 

3 Our convention for the coupling A differs from the one in Ref.^5] by a factor of 2. 
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which defines a self-consistency condition or gap equation. 
For Aij and m^ one uses the O(N) invariant ansatze 
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(2.5) 
(2.6) 



so that the equations for the effective masses separate as 



30 2 - w 2 + 3A a + (iV — 1)A, 



-i; 2 + A a + (iV + l)A 7r 



The gap equations become 
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and the effective potential takes the form 
V eS (m 2 a ,m 



4,< 



A 
I' 



2 ) 2 + V— (m 2 ,m 2 ,0) 



r2PPI 



A 



~ (3A 2 + (N 2 - 1)A 2 + 2(iV - 1) A CT A 7 



(2.7) 



(2J 



(2.9) 



where V 2PPl is the "quantum" part of the 2PPI effective potential. As has been shown in 
Refs. |1U[ ITS] these equations can be properly renormalized and the renormalized equa- 
tions have the same form. We do not discuss this here. For the numerical calculation 
we have used the renormalized versions of these equations; we have not put renormaliza- 
tion conditions but used an MS prescription. The renormalization scale fx refers to this 
prescription. 



3 Computation of the effective potential 

The basic relations given in the previous section can be used to compute the effective 
potential. We would have to solve the coupled system of gap equations and to insert the 
result into the 1PI effective action. This would imply that we would not only have to 
evaluate the sunset graphs, but also their derivatives with respect to m 2 and m\ . Here we 
prefer to work with an effective potential that leads to the gap equations by finding the 
extremum (maximum) with respect to variations of m 2 and m 2 . Instead of solving the gap 
equations whose algebraic and analytic form is already quite involved, we then can simply 
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use numerical algorithms for extremizing a function of two variables j2Zj- To this end we 
solve Eqs. ()2.7|) with respect to A CT and A n and insert these expressions into Eq. ()2.9j) . 
We denote this new potential by U(m 2 ,ml,4>). It can easily be verified, that Eqs. ()2.7|) 
again follow by extremizing this potential with respect to m\ and m 2 . The 1PI effective 
potential SIS db function of (p alone is obtained as 

V£\cf>) = U{ml,ml,<i)) (3.1) 

where m a and m, are the values which extremize (maximize) U for a given <f>. This 
procedure, as introduced by Nemoto et al. [B] in the Hartree approximation, generalizes to 
the case where higher order contributions are included into r^ PPI (see also Appendix 1X1 2). 
Here we include the two-loop contribution, the sunset diagram, as has been done previously 
for the N — 1 model by Smet et al. |16j . 

With these preliminaries we can now give our explicit equations: We decompose the 
potential U{m 2 r , m 2 , 0) into three parts: 

U = f/ c l ass + Z7i —loop + ^sunset • (3-2) 

The classical potential has the form (see jH]) 

^ci ass = \ml<f> 2 - ^ - 2A 1 + v 2 [ml + (N - l)ml) (3.3) 

~ %x{ N + 2) [ (iV + 1)m ' + 3(N ~ 1)m " " 2(N ~ 1)m ° ml + 2NX2V '- ' 

one easily checks that it takes its maximum if m 2 = \(3<f) 2 — v 2 ) and m\ = A(</> 2 — v 2 ) . The 
1-loop part is given by the "In det" contributions. At finite temperature these include the 
free energies, so the one-loop part of the effective action reads 

+Tj In [1 - ex P (-S CT (k)/r)] (3.4) 



+(iV -1)tJ-^L In [1 - exp(-^(k)/T)] 



(27T) 

In computing the sunset diagram it is convenient to decompose the finite temperature 
propagators into a zero temperature and a finite temperature (thermal) part proportional 
to 5(k 2 — m 2 )/ exp(— E^/Zc^T). The contribution of the sunset diagrams then consists of 
three parts (see e.g. Ref. 38 ) 

77 _ tt(°) i ttW i rr( 2 ) (o z\ 

'-'sunset w sunset ' u sunset ' w sunset ) 

with the T = contribution U^ nBet , the diagrams with one thermal line f/ S uLet an d the 

(2) 

diagrams with two thermal lines £/ sunse t, see Fig. 1. The T = part is given by 
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(a) Zero temperature part. 




(b) Graphs with one thermal line. 




(c) Graphs with two thermal lines. 

Figure 1: Contributions to the sunset diagram at finite temperature. Solid lines represent 
the zero temperature parts of the propagators, finite temperature parts are denoted by 
interrupted lines. 
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u. 



(0) 



[3I aaa + (N- l)I a ^\ 



(3.6) 



it is represented graphically in Fig. The diagrams with one thermal line, see Fig. ^p, 
contribute 



U. 



(i) 

sunset 



-A 2 2 \9lL a + (N-1) (2lf , + l£. 



(3.7) 



The symbols for the thermal lines are underlined. Similarly the diagrams with two thermal 
lines, see Fig. Ht, contribute 



(2) 



sunset 



X 2 (b 2 \9I p , + (N — 1) (2/ /3 , +J /3 I 



(3.8) 



The precise definition of the Feynman integrals 1^ , J^-a and 7^ fc with zero, one and two 
thermal lines, respectively, as well as their analytic form are presented in the Appendices. 
It is understood that their divergent parts are removed. 



4 Discussion of the numerical results 

As we have stated previously we do not solve the two coupled gap equations but instead 
we maximize the potential U{m 2 T) m 2 , 0). We present our numerical results for the case 
iV = 4 with A = 1 and A = 0.1 4 . The mass scale is fixed by taking v = 1, and we choose 
the MS renormalization scale ~p = 1. In Fig. El we display the value of (f) at the minimum 
of the effective potential, the thermal expectation value which we denote by v(T). If we 




Figure 2: The expectation value v(T) of for iV = 4, \i = v = 1, A = 1 (squares) and 
A = 0.1 (diamonds). 

4 The choice A = 1 corresponds to A = 2 in the normalization of Ref. ^U] an d to A = 6 when compared 
to the N = 1 study of Ref. [El 
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choose A = 1 we see a phase transition towards the symmetric phase v(T) = for T > T 
with T ~ 1.7. For A = 0.1 the critical temperature is about T ~ 1.475. 

The behavior of the effective potential as a function of <p near the phase transition is 
displayed in Fig. 01 which clearly indicates that the phase transition is of second order. The 
behavior, for the same parameters (A = 1) but without the sunset diagram (i.e. Hartree 
approximation, see Ref. [Tfi]). is displayed in Fig. HJ The two minima characteristic of a 
first order phase transition are well visible. It is well known that a phase transition of first 
order is found in the Hartree approximation. As apparent from the scale on the y axes 
and from the tiny temperature range, Figs. 3 and 4 represent "microscopic" pictures of 
the two phase transitions. 

The temperature dependence of the sigma mass M a as defined by the curvature of the 
effective potential at its minimum is shown in Fig. As to be expected it goes to zero at 
the phase transition temperature, the zero is approached linearly if one plots Mj. 

In Fig. El we display the temperature dependence of the variational masses m a and 
m-n at the minimum of the effective potential, <ft — V {T)- The variational sigma mass mo- 
behaves similarly as the sigma mass M a obtained from the effective potential. The mass 

becomes identical to the mass m a above the phase transition, but does not vanish below 
the phase transition. It was found already in the one-loop analysis that, as a "violation 
of the Nambu-Goldstone theorem", the self-consistent pion masses do not vanish when 
the symmetry is broken. It has been argued that these self-consistent masses are not the 
physical pion masses; indeed they are not: they are variational parameters and the effective 
potential Kff(0) has of course N — 1 flat directions at its minimum. But the hope - or 
expectation - that the discrepancy between the pion mass as computed from the effective 
potential, and the pion mass as a variational parameter would disappear, turns out to be 
fallacious. Indeed there is a simple physical reason for this result, as will be discussed 
below the next paragraph. 

As the thermal integrals require both masses to be real we have looked for the extremum 
with respect to m a and m n . So in our numerical approach neither w? a nor can get 
negative and the well-known instability which occurs for AT = 1 in the region where the 
potential has a negative curvature is avoided by fiat, the maximum simply occurs at the 
boundary of the "physical region" of real pion and sigma masses. It has to be said, though, 
that in this case we do not solve the gap equation which in fact becomes meaningless. In 
the large-A^ limit such a construction leads to an effective potential that is flat in the 
region <fi < v , and therefore convex. We do not want to enter into this discussion here, 
we simply state that these regions require another approach and that we have to discard 
them. Such regions occur at low temperatures only, and of course they do not include the 
region around the minimum of the effective potential. At higher temperatures, but well 
below the phase transition the effective potential has regions of negative curvature, but 
both m a and m n , and therefore the effective potential are still real, as the variational mass 
m a is not equal to the curvature of the potential. The parameter m n , which is imaginary 
below the minimum in the large- A^ approximation, also becomes real for all values of <fi at 
temperatures well below the phase transition. 

However, we are faced with an even more important new feature: the fact that the sigma 
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(a) A = 1 for temperatures T = 1.62, 1.66, 1.69, 1.70 and 1.72 
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Figure 4: The behavior of near the critical temperature in the Hartree approxi- 

mation; the curves are for T = 1.46, 1.47, 1.48 and 1.49; parameter set as in Fig. Elwith 
A = l. 




Figure 5: The sigma mass M a obtained from the effective potential as a function of tem- 
perature; parameter sets as for Fig. El 
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Figure 6: The variational masses m a (diamonds) and m n (squares) as functions of temper- 
ature; parameter sets as for Fig. |21 



can decay into two pions if m a > 2m^. The sunset diagram with one thermal sigma line 
and two pion lines acquires an imaginary part in this case. In our computations we have 
simply omitted this imaginary part, but obviously we would not be able to solve the gap 
equation in regions where such a decay is possible. While we did not exclude these regions 
we have to consider the real part of the effective potential in these regions with suspicion. 
In contrast to the problem of imaginary masses and the associated instability these regions 
do include the minimum of the effective potential for temperatures up to almost the phase 
transition temperature. We display the (trial) masses m a and m n at the minimum of the 
effective potential in Fig. El Around the phase transition itself we find m a < 2m n , so 
that the behavior of the effective potential in the critical region, as plotted in Fig. El is 
not affected, but the results below T ~ 1.5 for A = 1 and below T ~ 1.4 for A = 0.1, 
respectively, have to be taken with some caveat. 

This finding has important consequences: Of course an unstable particle can coexist 
with its decay products at finite temperature, but this situation requires an approach where 
the transitions are taken into account; however, this is not the case in this approximation, 
and indeed with the entire formalism used here. Indeed in the regions affected by this 
instability our approximation becomes inconsistent, and this should be so a fortiori if one 
considers the massless physical pions. As the masslessness of the Goldstone particles is 
an important aspect of spontaneously broken symmetry, this problem should be studied 
in detail. Of course in the applications to real pions in the linear sigma model the pions 
receive a finite mass due to explicit symmetry breaking, and the sigma particle is considered 
usually as being of a problematic status anyway, hinting at the limitations of the model as 
an effective theory of strong interactions. 
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5 Summary and Outlook 



We have analyzed here the 0(N) linear sigma model in the 2PPI formalism beyond the 
leading order, in which it coincides with the Hartree approximation. As in the Hartree 
approximation and in the N = 1 version we find that the effective mass of the pion 
quantum fluctuations is different from zero in the broken symmetry phase, so that a naive 
particle interpretation, suggested by the large-iV analysis, becomes problematic. As in 
the N — 1 version of the model ^Sj the phase transition, which is first order in the 
Hartree approximation, becomes second order. In addition to the N — 1 case there is 
a new instability associated with the possibility of the decay a — > 2tt. This will not be 
problematic at low temperatures and for small couplings, but whenever the sunset diagrams 
become important it requires reconsidering the entire framework. We find that near the 
phase transition the sigma fluctuations become stable, as they are trivially in the symmetric 
phase. 

Our analysis should have some bearing on nonequilibrium simulations as well. The 
non-vanishing effective mass of the "Goldstone" quantum fluctuations makes it hard to 
maintain a naive particle interpretation; but this is the case a fortiori for any nonequi- 
librium simulations that include higher order diagrams, for approximations in which the 
propagator is not an effective free particle propagator. In addition, however, it becomes 
obvious that the additional instability that occurs only for N > 1 will lead to other and 
new aspects of such nonequilibrium simulations, when compared to those for the large- N 
case. While it is certainly important to understand thermalization, the instabilities both of 
the 1-loop effective potential as those introduced by particle decay may have consequences 
of similar importance, and conclusions drawn from N = 1 simulations may therefore miss 
essential aspects for models with spontaneous symmetry breaking. 

We would finally like to remark that all the complications found in thermal equilibrium 
occur in nonequilibrium studies as well, both in the preparation of the initial state and in 
the analysis of the final state, as well as in renormalization. Therefore it is mandatory that 
such equilibrium studies are being pursued in parallel to the nonequilibrium ones. 
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A Comparison of 2PPI and 2PI effective action 



We compare the 2PI with the 2PPI effective action using a loop expansion. The Lagrangian 
of the classical action is given by Eq. (|2.1jl but here we set N = 1 for the sake of convenience. 



A.l 2PI expansion 

The 2PI effective action reads (cf. Eq. (2.9a) of Ref. [TB] ) 

T 2Pl [<J),G} = S[<J)} + i^lndetG- 1 + i^Tr V^G + T 2 2 PI [<f),G} , (A.l) 

where T^ 1 contains all higher 2PI corrections; the first graphs of the loop expansion are 
shown in Fig. 




Figure 7: Two-loop and some (not all) three-loop contributions to the 2PI effective action: 
double-bubble, sunset, basketball and another three- loop graph. A line represents the full 
propagator as defined in Eq. ([A. 21) . 



The propagator G takes the form 

i 

where the equation for the proper self-energy £ follows from the condition 5T ^ G = 

So the 2PI self-energy is obtained by cutting a line of a 2PI vacuum graph and considering 
a combinatorial factor. Again, the leading terms of the self-energy in a loop expansion are 
displayed in Fig. |H1 



A. 2 2PPI expansion 

The 2PPI effective action, proposed by Verschelde and Coppens [Tl], is a variant of the 
"effective action of composite operators" by Cornwall, Jackiw and Tomboulis ^Hj. We 
briefly repeat the formal derivation in Ref. ^1] without going too much into detail. 

In the 2PI formalism one deals with a bilocal composite operator $(x) $(y) which is 
coupled to an external (bilocal) source K(x,y), while in the 2PPI approach one keeps this 
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Figure 8: Some one- and two- loop contributions to the proper self-energy of the 2PI expan- 
sion: tadpole, fish, sunset and another two-loop graph. A line represents a full propagator 
as defined in Eq. (|A.2|) . 



source local by construction. Using Euclidean space-time, as in Ref. [14;, one defines the 
following effective action of local composite operators 

T 2PPI [0,A] = W[J U J 2 ]- J!-0-^J 2 -(0 2 + A) (A.4) 



where 



5W 



and 



,)./, ($) = (A - 5) 

exp{-W[Ji, J 2 ]} = jv$ exp{-(5 , [$] + Ji • $ + J 2 • $ 2 )} . (A.7) 

The external fields (sources) J\ and J 2 are both local and fix the expectation value of $ 
and <3> 2 . The effective equations of motion turn out to be 

SMI = (A.8) 

0(p 

<5A " "2 Ja - (A ' 9) 

We do not want to explain in detail the combinatorial trick used in Ref. ^4] for the 
derivative — — S J;^' A ^ in order to sum all 2PR graphs. Eventually, the result is the complete 
effective action of the 2PPI formalism. It consists of the classical action, the "quantum" 
part and a constant which prevents double-counting: 

r 2PPI [0, A] = S[<j>] + r 2PP V, M 2 ] - ^3 AA 2 . (A.10) 

The "quantum" part r^ PPI Q f the 2PPI action contains the one-loop "In det" contribu- 
tion and all 2PPI graphs without the "double-bubble" (see Fig. EJ). A two-particle-point- 
irreducible (2PPI) graph is 1PI and stays connected whenever two internal lines meeting at 
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rf PI = ( ) t 1 f— ^ - i-t-T 

Figure 9: "Quantum" part of the 2PPI effective action including some three- loop contri- 
butions. Here, a line represents a bubble-resummed, i.e., local, propagator as defined in 
Eq. ([A. II)) . Note that the last graph does not appear in T 2 , pl because it is 2PPI but 2PR. 



the same point (vertex) are cut [Hj . These graphs are to be computed using a propagator 
with an effective mass Ai 

G « = ■ < A ' n > 

Since this propagator always remains local - even if the exact 2PPI effective action is 
computed - it will never be equal to the physical two-point Green function. 

The effective mass consists of the "classical" mass (-Aw 2 ), the seagulls \<j) 2 and the 
local self-energy A 

M 2 = \(3(j) 2 - v 2 ) + 3 \A . (A.12) 

The self-energy is — like in 2PI — obtained as a derivative of the "quantum" part of the 
effective action 

This is equivalent to cutting a line (by deriving with respect to the propagator G(k)) 
and then connecting the two endpoints to a common third point by considering the inner 
derivative 

6 r 5G(k) 5 r -1 5 

SM 2 ' ~ J m 5M 2 5G(k) ~ J [ 1 (k 2 + M 2 ) 2 5G(k) ' { } 

The first (one-loop) "quantum correction" to the 2PPI effective action only consists of the 
"In det" term in T 2PPl . This approximation is equivalent to what is called "Hartree" in 
2PI, namely the resummation of tadpoles or daisy and super-daisy graphs. The first two- 
loop contribution to the 2PPI effective action is the sunset diagram. The mass corrections 
resulting from all one-, two- and three- loop vacuum diagrams of Fig.lHlare shown in Fig.fTUl 



A = + 





Figure 10: Contributions to the 2PPI self-energy with some three-loop contributions. A 
line represents a local bubble-resummed, propagator as defined in Eq. (jA.llJ) . Note that 
the last two graphs do not appear in T 2PI because they are 2PPI but 2PR. 
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Now, we use Eq. (jA.12|) to express A in terms of M 2 and insert this expression into the 
effective action (jA.lOjl . We finally obtain the 2PPI effective action in terms of the mean 
field and the effective mass M.. Restricting to a homogeneous mean field and using the 
explicit form of the classical Lagrangian (]2.1j) the effective potential reads 

V cS (<p,M 2 ) = \m 2 4> 2 -^ 

+V 2PP \<P,M 2 ) -^- X [M 2 + Xvf . (A.15) 

This is the effective potential from which the equations of motion can be (re)-obtained by 
varying (f) and M 2 . This procedure can be generalized to obtain an effective action for 
nonequilibrium dynamics and a conserved energy functional [39J. 

B The sunset diagram at T = 

The sunset diagram at T = has been evaluated previously by several authors j47H l4"T] . 
Here we are interested in this diagram with external momentum zero, and with at least 
two internal lines of equal mass. We give here some technical details. 

B.l The a a a diagram 

The sunset diagram for equal masses has been given in [40] . The authors use d = 4 + e 
and omit the factor l/(27r) d . Alternatively one may use the expression given in jH] for 
different masses and set all masses equal. These authors use d = 4 — 2e and likewise omit 
the factor l/(2ir) d . Accomodating both expressions to our standard 

r d A ~ £ p r d*- £ g 1 

aac J (27r) 4 ^7 (27r) 4 - e (p 2 + M 2 )(g 2 + M 2 )((p-g) 2 + M 2 ) 1 ' ' 



we find 



.2 2 / 4vr/i 2 

laaa = 71 77 W + ~ ln ~ 777T ~ IE + ~ 




where C1(0) is the Clausen function 



CW = (B.3) 



k=l 
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B.2 The airir diagram 

Defining the sunset integral for one a line with mass M = m a and two pion lines with 

mass m = m, as 

d 4 ~ £ p f d 4: ~ e q 



1 



(2vr) 4 " e J (2vr) 4 - e (p 2 + M 2 )(q 2 + m 2 )((p - q) 2 + m 2 ) 



2\ ' 



(B.4) 



we again use the expressions given by Davydychev and Tausk jHj which we have cross- 
checked with the the reduction formula of Ref . 02] . When adapted to our conventions we 
obtain 

2m 2 (l + 2z)[2 2/ 4tt// 2 3 2z , \ 

J - = \7 2 + ir^ 2 - ~ lE + 2~ TTT z lnAz 



2 2 5 



/ 47T/i 2 3 2Z , . \ 7T 

+ In — V-7i?H In 4,2 H h- 

l M 2 2 1 + 2z y 12 4 



^ 1 -^). ln 2 4z _V^ 1 -^ 



1 + 2^ 



l + 2z 



C%) 



where 2 = M 2 /4m 2 and 9? = arccos(l — 2z) 



(B.5) 



C The sunset diagram with one thermal line 

If one of the lines of the sunset diagram is replaced by a thermal line it takes the form 

d 3 ~ e p g , f d A ~ e q 1 



I 13 



(27r) 4 - {q 2 + m 2 ){{p - q) 2 + m 2 ) 



(27r)3-^ fc (p) 

Here E k (p) = J p 2 + m| and (p) is the Bose- Einstein distribution function 



(C.l) 



4(p) 



1 



exp(E k /T) - 1 • 



(C.2) 



The second integral is the fish diagram with the external Euclidean momentum p which is 

^2 _ ,^2 



on shell, i.e., p = —m k 

F(mi,mj;m k ) 
Explicitly one finds 

F(m h m j ]m k ) ■■ 



d 4 ~ e q 



(2vr) 4 - e {q 2 + m 2 ){{p - q) 2 + m 2 ) 



(C.3) 



mi 



(4tt) 2 [e 

— / da\n 
Jo 



- -7 B + ln47r-ln— ^ 



777, 2 777 2 

a — \ + (1 — a) — ij- — a(l — a) 

m k m k 



(C.4) 
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The a integration can be done analytically, leading to various expressions in terms of 
logarithms or inverse trigonometric functions, depending on the relations among the three 
masses. We note in particular that for m k > + rrij the integral (|C.4J) develops an 
imaginary part. In the present context this happens for m k = m a and = rrij = m, < 
m a /2. The thermal integral 7^., fe thereby gets an imaginary part as well, reflecting the fact 
that in the heat bath the sigma particles can decay into or be produced by pions. The 
sunset integral with one thermal line factorizes: 

l tj\k = F ( m i> m f> m k )I^_ e {m k ) (C.5) 

with 

/ -M=/(syQbr* (p) ' (a6) 

The divergent part obviously takes the form 

'm = (4^- t M • (C.7) 

The general expression agrees with Ref. [SB], Eq. (3.7), we use a different regularization, 
however. 



D The sunset diagram with two thermal lines 

We define the sunset diagram with two thermal lines as 



i\jk 



d 3 ~ 



P 



(2vr)3-2E fe (p 



-"fc(p) 



d 3 ~ 



(27r)3-2^(q; 



r,s=± 



(Pr + Qs 



■HI; 



(D.1) 



Here p r and g s are Euclidean momenta with p± = ± (z£?fe (p) , p) and q± = ± {iEj (q) , q) . The 
integration over the angle between p and q can be done analytically, with the result 



i\jk 



32tt 4 Jo E k (p 



pdp p 
<(P) 



qdq p 

nAq) In 



o Ej(q 



(D.2) 



where 



Y ± = {[E k (p) + EM)? ~ Ef(p ± q)} {[E k (p) - EM)? ~ Ef(p ± ?)} ■ (D.3) 



The integrand has logarithmic singularities within the region of integration which have to 
be treated with care in the numerical integration over p = |p| and q — |q|. 
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